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Gravitational microlensing as a test of stellar model atmospheres 
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Abstract. We present calculations illustrating the potential of gravitational microlensing to discriminate between 
classical models of stellar surface brightness profiles and the recently computed "Next Generation" models of 
Hauschildt et al. These spherically-symmetric models include a much improved treatment of molecular lines in 
the outer atmospheres of cool giants - stars which are very typical sources in Galactic bulge microlensing events. 
We show that the microlensing signatures of intensively monitored point and fold caustic crossing events are 
readily able to distinguish between NextGen and the classical models, provided a photometric accuracy of 0.01 
magnitudes is reached. This accuracy is now routinely achieved by alert networks, and hence current observations 
can discriminate between such model atmospheres, providing a unique insight on stellar photospheres. 
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1. Introduction 

Recently gravitational microlensing has been shown to 
be an effective tool for stellar astrophysics, primarily 
through high time resolution observations of extended 
source events - i.e. where the lensed star cannot reason- 
ably be approximated as a point source. 

The microlensing signatures of extended sources were 
initially discussed in Bontz (1979), Gould (1994), Witt & 
Mao (1994) and Nemiroff & Wickramasinghe (1994) for 
the case of a star of uniform surface brightness, where 
the finite size of the star produces a significant deviation 
from the point source light curve. The primary motivation 
for the authors in the 1990's was to use extended source 
effects to better constrain the properties of the lens - thus 
distinguishing between e.g. LMC and galactic lenses. 

The microlensing signatures of a non-uniform disc 
was first considered by Valls-Gabaud (1994), Bogdanov & 
Cherepashchuk (1995), Witt (1995) and Simmons, Willis 
and Newsam (1995), who showed that the lightcurves 
would display a chromatic signature as the lens effectively 
sees a star of different radius at different wavelengths. 
Valls-Gabaud (1994, 1998) modelled linear, quadratic and 
logarithmic limb darkening, with coefficients computed for 
the Johnson U to K wavebands using ATLAS LTE model 
atmospheres from Kurucz (1994), and also predicted a 
spectroscopic effect, where the line profiles would change 
during the microlensing event, providing another probe of 
atmospheric structure. This was confirmed by Heyrovsky, 
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Sasselov & Loeb (2000) for a particular red giant model 
atmosphere. 

Gould (2001) reviews recent observational progress 
in detecting non- uniform surface brightness effects. The 
clearest evidence of limb darkening to date comes from 
analysis of event MACHO 97-BLG-28 (Albrow et al. 
1999). This was an example of a binary lens event with a 
cusp crossing - arguably the most favourable type of event 
for stellar astrophysics studies. The photometry and sam- 
pling coverage for this event were of sufficient quality to 
fit a 2-parameter limb darkening law in V and I. 

Although observations already strongly favour limb 
darkened atmospheres over uniform brightness models, 
the use of microlensing as a discriminant between limb 
darkened models is a relatively new subject. Albrow et 
al. (2001) present a detailed analysis of event OGLE 99- 
BUL-23, comparing the V and I light curves with several 
different model atmospheres, including for the first time 
a treatment of the errors in the lens parameters and their 
correlation with the estimated limb darkening coefficients. 
Notwithstanding the rigorous statistical approach of that 
paper, which 'raises the benchmark' for future analyses, 
only linear limb darkening models were considered. The 
current 'state of the art' atmosphere models are the 'Next 
Generation' models of Hauschildt et al. (1999a, b), which 
include a much more detailed treatment of molecular ab- 
sorption lines in the atmospheres of cool giants and con- 
sider spherically symmetric atmospheres, as opposed to 
plane parallel models used previously. The aim of this 
Letter is to compute the microlensing signatures of typical 
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extended source events, lensed by point and fold caustics. 
This allows a comparison of the differing stellar atmo- 
sphere models to be made and so, stimulating interest in 
observations of extended source events as a test of model 
atmospheres. 

2. PHOENIX Next Generation model 
atmospheres 

Computations of stellar surface brightness profiles have 
been carried out for several decades but until very recently 
were based on approximate fits to the limb darkening, that 
is the dependence of the intensity, /(/it), as a function of 
(the cosine of) emergent angle, fi. For instance a simple 
linear model, namely 7(/z) — Io [1 — ci(l — /x)] was often 
sufficient, while the detailed study of eclipsing binaries 
sometimes required two parameter models, such as the 
"square root" /(/j) = Io [l — c 2 (l — /i) — c 3 (l — y/JI)] or 
"logarithmic" model 7(/z) = Io[l — c 4 (l — /i) — cs/xln/x]. 
Each coefficient Cj depends on the temperature, gravity 
and chemical composition of the source, but the range of 
validity of these "laws" is often very limited (see Valls- 
Gabaud 1998 and references therein). For instance the 
square root formulation provides a good fit for hot stars, 
while a quadratic equation would be enough for cooler 
stars. A new non-linear formulation, using 4 parameters 
and valid across the entire range of effective tempera- 
tures and surface gravities has recently been suggested 
by Claret (2000). Claret (2000) also made non- linear fits 
to the ATLAS model atmospheres, improving upon the 
previous calculations made by Van Hamme (1993). 

By contrast, the recent PHOENIX "Next Generation" 
stellar atmosphere models, computed by Hauschildt et 
al. (1999a, b), considerably improve upon these models in 
several important respects, and in particular by carrying 
out calculations assuming spherical geometry. Moreover, 
the intensity calculations are based on a huge library of 
atomic and molecular lines, with about 2 x 10 8 molecular 
lines contributing to a typical giant atmosphere model at 
T cff = 3000 K. Claret (2000) also fitted his new non-linear 
limb darkening equation to some PHOENIX models. 

The dramatic difference in the dependence of limb 
darkening on emergent angle between the ATLAS models 
and PHOENIX models is illustrated in Figure 1, which 
shows the intensity profiles for a giant star of T e g = 4250 
K and log g = 0.5, in four Johnson bands: V, R, I and K. 
The solid curve shows the PHOENIX profiles, while the 
dashed, dash-dotted and dotted curves denote the linear, 
logarithmic and square root models respectively. It is im- 
mediately clear that there is a dramatic decrease in the 
predicted intensity of the PHOENIX models as one ap- 
proaches the limb of the star, i.e. at \i ~ 0.2. This feature 
does not appear in ATLAS models and arises from the 
effects of spherical symmetry, which includes finite opti- 
cal depth for rays close to the limb - a treatment which 
is absent from ATLAS models which utilise plane parallel 
atmospheres in which the optical depth of limb rays is still 
finite and so provide significant intensity out to fi = 0. The 






Fig. 1. Intensity profiles for a giant star of T e e = 4250 K and 
log g = 0.5, in four Johnson passbands: V, R, I and K. The 
solid curve shows the PHOENIX profiles, while the dashed, 
dash-dotted and dotted curves denote the ATLAS linear, log- 
arithmic and square root models respectively. 

question then arises: is microlensing sufficiently sensitive 
to detect this limb feature in the atmospheres of extended 
sources, and thus to test the NextGen models against real 
observations? 



3. Computing microlensing lightcurves 

For a point source the amplification due to a point mass 
lens is A = (u 2 + 2) / (uy/u 2 + 4) , where u(t) is the im- 
pact parameter, the projected angular separation between 
lens and source, measured in units of 9e, the angular 
Einstein radius of the lens, defined as 
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where Di and D s are the lens-observer and source-observer 
separations respectively and M the mass of the lens. If the 
impact parameter is comparable to the angular radius of 
the source, the point source amplification function breaks 
down and it becomes necessary to calculate the amplifica- 
tion as an integral over the source star, given by 



A v {t) = 



J J I„(r, 9) A(r,6,t) rdrd9 
J J I v {r, 9) rdrdO 



(2) 



Thus, the microlensing lightcurve contains unique infor- 
mation on the source surface brightness profile, I v (r,9), 
which is a far more sensitive test of the model atmosphere 
than its integrated value, the emergent flux. 

In the case of binary lenses, the amplification func- 
tion takes a more complex form, although for an extended 
source it is again given by an integral over the source. In 
this Letter we are only concerned with so-called fold caus- 
tic crossings, and so we adopt the 'square-root' approxi- 
mation for the amplification inside the caustic structure 
(Schneider, Ehlers and Falco, 1992) which is valid within 
a few angular source radii of the caustic and takes the 
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Fig. 2. A comparison of the V, R, I and K lightcurves pro- 
duced by the transit of a point lens with minimum impact pa- 
rameter uo = 0.0 of a 4250 K, log g = 0.5, 0.1 6e source, mod- 
elled with PHOENIX and linear limb darkening. Magnitude 
differences are computed according to Equation w in the text. 
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Fig. 3. A comparison of the V, R, I and K lightcurves as in 
Figure ^, but now for an event with uo = 0.09. 

form A(x) = Ao + l/y/x. Here Ao is the total amplifica- 
tion of other images outside the caustic, which we assume 
to be constant during the caustic crossing, x is the per- 
pendicular distance in units of the angular Einstein radius 
of the binary lens, from an element of the source to the 
caustic. Note that for elements of the source outside the 
caustic structure there are no additional images formed, 
so A(x) = A . 

4. Results 

Figure [2] shows a comparison between two sets of mi- 
crolensing lightcurves for a point lens transit event: one 
set calculated using a linear limb darkening law and pa- 
rameters from the ATLAS grid by Kurucz (1994) (see 
Valls-Gabaud, 1998); and the second set calculated us- 
ing the PHOENIX Next Generation models as discussed 
in Section 2. In this example we compare the models 
for a source with of angular radius 0A9e, logg = 0.5, 
T e ff = 4250 K and solar metallicity. The abscissa shows 
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Fig. 4. A comparison of the V , R, I and K lightcurves pro- 
duced by fold caustic crossings of 4250 K, logg = 0.5, 0.1 6e 
sources with PHOENIX and linear limb darkening according 
to Equation |j| 

the time evolution of the event, in units of the Einstein 
crossing time, ts, of the lens - i.e. the time taken for the 
lens to move 1 angular Einstein radius, to denotes the time 
at which the impact parameter is a minimum, uq] in this 
example we take uq = - i.e. the lens transits through the 
centre of the source. The ordinate shows the magnitude 
difference between these models, defined as 

magncxtgon-maglincar = 2 -5 log ffT^J ~ 2 - 5 lo § \^~) ^ 

where Fll , -Ful > -Fln and Ftjn are the lensed and unlenscd 
linear limb darkened, and lensed and unlenscd NextGen 
fluxes respectively. Four colour bands are used to illustrate 
the strong chromatic differences between the models: V, 
R, I and K bands, represented by straight, dashed, dash- 
dotted and dotted lines respectively. 

Two main features are immediately apparent in Figure 
g. Firstly, just before and after the lens transits the 
source there are positive "spikes"; these occur as the 
PHOENIX/NextGen model is considerably more limb 
darkened than the linear model close to the limb. However 
these spikes are very narrow. The second feature is the 
broad 'brightening' during the central phase of the transit; 
this also occurs as a consequence of the strong NextGen 
limb darkening, which makes the source appear smaller as 
the lens crosses its centre, causing a larger amplification. 
Thus, if a source with a NextGen type atmosphere is fitted 
by a linear limb darkening model, the source size would 
be systematically underestimated. 

Figure || compares lightcurves for the same source as 
Figure ^ but with an impact parameter of 0.09 9e, i-e. the 
lens 'just' transits the source. In this case the only features 
are the two upward spikes, due to the NextGen intensity 
being considerably smaller than the linear model intensity 
at the limb of the star. Again, however, these spikes are 
very narrow and around minimum impact parameter the 
difference between the models is much smaller, since the 
amplified flux is dominated by a region slightly closer to 
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the centre of the source where there is less difference be- 
tween the models (see Figure [I] at /j = 0.9). In particular, 
no brightening effect is seen since the lens does not cross 
the central region of the source. The amplitude of the ef- 
fect is, however, similar: at least 0.01 mag, which is easily 
detectable. 

Figure |] compares the NextGen and linear models dur- 
ing a fold caustic crossing event, shown from when the 
source is 3 source radii outside to 3 source radii inside 
the caustic structure. Here the source has T e g = 4250 
K, logg = 0.5 and solar metallicity, but now has an an- 
gular radius of 0.01 0e - a more realistic value for e.g. 
typical Bulge events, which results in larger amplification 
changes. The most striking feature on Figure f| is the large 
spike as the source begins to enter the caustic; here the 
flux is dominated by the small region of photosphere very 
close to the limb and directly beneath the caustic - the 
surface brightness of which varies dramatically between 
the two models. As the event proceeds the magnitude dif- 
ference becomes much smaller, particularly in mid-transit 
when the limb differences between the models are diluted 
by the larger intensity from the central part of the source. 
However, we again see a broad 'brightening' effect as in 
Figure ||, due to the extreme NextGen limb darkening pro- 
ducing an apparently smaller source. As the trailing limb 
crosses the caustic we see only a small peak since there 
is already considerable amplification from the rest of the 
photosphere. In practice, of course, we would expect to ob- 
serve these events in reverse: it is easier to predict when 
a source will exit a caustic than when it will enter! 

5. Conclusions 

Figures 2-4 clearly suggest that it is within the capabili- 
ties of current intensive microlensing monitoring programs 
to discriminate between the chosen atmosphere models 
from a well parameterised lightcurve. For the case of a 
point caustic event, with e.g. timescale tg — 15 days, the 
transits shown in Figure 2 and 3 would last for about 3 
days, with the sharp features at each limb crossing a few 
hours in duration. The sampling rates regularly achieved 
by e.g. the PLANET collaboration during recent events 
would easily suffice to resolve these limb crossing features 
with a reasonable number of data points. It must be noted, 
of course, that since Figures 2 and 3 are for point lens 
transits, the probability of such an event - and indeed of 
its immediate detection - is small. Given that the event 
is alerted during the rise phase as a potential extended 
source, however, the detection of the limb crossing spikes 
would be highly likely. 

The situation for the fold caustic event is more en- 
couraging. Again, the timescale for the caustic crossing is 
of the order of several hours, which is typical of "lensing 
anomaly" events already regularly monitored. Moreover, 
as pointed out previously, the chief advantage of fold caus- 
tic crossing events is the predictability of the source ex- 
iting the caustic structure, thus allowing intensive mon- 
itoring to be scheduled in advance. As the trailing limb 



of the star leaves the interior of the caustic one will see a 
dramatic difference between NextGen and the other mod- 
els since it is effectively only the (strongly darkened) limb 
of the star which is being amplified. The signature of this 
limb crossing is seen in Figure 4 to be of order 0.2 mag. 
in Johnson B, V, R and I - which would render it very 
easily detectable with existing photometric precision. 

To compare any stellar atmosphere models using real 
data would involve simultaneously determining best-fit 
lens parameters for the event under study. As remarked 
previously, Albrow et al (2001) have included the effects of 
lens parameter errors in estimating limb darkening coeffi- 
cients. We will carry out a similar numerical study, using 
realistic simulated observations, in future work and at a 
more advanced stage one must also consider the possibil- 
ity of spots further effecting the lightcurve (see Hendry 
et al. 2002). It seems clear, however, that the magnitude 
of the difference between NextGen and ATLAS models 
is eminently detectable - even allowing for uncertainties 
in the fitted lens parameters. This is an essential step in 
the determination of fundamental stellar parameters. For 
instance, high accuracy measures of stellar radii via eclips- 
ing binary light curves rely critically on these model at- 
mospheres (e.g. Orosz and Hauschildt 2000) as do inter- 
fcromctric measures (e.g. Davis et al. 2000, Wittkowski et 
al. 2001). Thanks to microlensing we have now almost di- 
rect access to the distribution of intensity across a source, 
and hence we can probe stellar photospheres with unprece- 
dented detail. 
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